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Abstract 

We study the thermodynamic behavior of branched polymers. We first study random 
walks in order to clarify the thermodynamic relation between the canonical ensem- 
ble and the grand canonical ensemble. We then show that correlation functions for 
branched polymers are given by those for </> 3 theory with a single mass insertion, not 
those for the </> 3 theory themselves. In particular, the two-point function behaves as 
1/p , not as 1/p 2 , in the scaling region. This behavior is consistent with the fact that 
the Hausdorff dimension of the branched polymer is four. 
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1 Introduction 

Branched polymers are the simplest generalization of the random walk and have been stud- 
ied extensively H p| . It is of great importance not only in statistical physics but also in 
particle physics, in particular, for understanding the critical behavior of random surfaces 
and quantum gravity@|||Q. Our recent interest in branched polymers arose in our attempt 
to formulate superstring theory nonperturbatively. In the paper |J, we have studied the 
dynamics of the type IIB matrix model in such a framework (JF| and see || for review). 
In our matrix model approach, the eigenvalues of matrices are interpreted as space-time 
coordinates. In these investigations, we find the system of branched polymers in a simple 
approximation. Although it is far from the flat four dimensional manifold, branched poly- 
mers share the same (fractal) dimensionality four with our space-time. It might be the first 
indication that superstring can explain the dimensionality of our space-time. 

In this paper, we comment on a field theoretic description of branched polymers. It is 
well-known that a system of random walks is described by a free scalar field theory if there 
is no effect of self-avoidance. Similarly it is widely believed that the system of branched 
polymers is described by a scalar field theory with a three-point coupling, that is, 3 scalar 
field theory. We will, however, show that it is not so by treating the universal part of the 
partition function carefully. The system of branched polymers without self-avoidance can 
be exactly solvable by introducing the grand canonical ensemble and using the so called 
Schwinger-Dyson technique. In order to extract the correct large N limit (N is the system 
size) or the thermodynamic limit, we have to check that the grand canonical ensemble is 
dominated by the larger size system. In other words, we have to extract the universal part. 

Our claim is that we need a single mass insertion in each m-point correlation function 
of the 3 scalar field theory in order to describe m-point correlation functions in branched 
polymers. Mass insertion here means a change of a propagator in each m-point function 
from an ordinary one, 1/p 2 , to 1/p 4 . In particular, the two-point function behaves as 1/p 4 , 
not 1/p 2 . Let us count the number of points which lie within distance R from a fixed point 
in d(> 4) dimensions. In the random walk, it can be estimated as R 2 by using the two-point 
function. We obtain R 4 for branched polymers by using the 1/p 4 type propagator. So our 
finding are consistent with the claim that branched polymers are four dimensional fractals. A 
multi-point correlation function is given by a sum of graphs of the corresponding correlation 
function for 3 scalar field theory at tree level with a single mass insertion in each graph. 

Our main results were announced in [ID|. In this paper we would like to give a fuller 



account of our results by providing more detailed derivations and explanations. The organi- 
zation of this paper is the following. In section 2, we review random walks in order to clarify 
the relation between the canonical and grand canonical ensemble. In section 3, we investigate 
branched polymers. First we define a canonical ensemble for a system of branched polymers 
(sec. 3.1) and then introduce grand canonical ensembles (sec. 3.2). We emphasize that the 
definition of a grand canonical ensemble is not unique. In sec. 3.3, we solve Schwinger- 
Dyson equations for the conventional grand canonical ensemble and obtain the results which 
correspond to the correlation functions of a scalar 3 theory. In sec. 3.4, we consider the 
thermodynamic limit of the correlation functions and obtain the correct universal behavior 
of them. In section 4, we give a physical interpretation why the propagator behaves 1/p 4 
in branched polymers. Section 5 is devoted to conclusions and discussions. We have two 



appendices A and B which derive the partition function and the two point function in the 
canonical ensemble of branched polymers. 

2 Random Walks 

In this section, we give a brief introduction to random walks in order to clarify the thermo- 
dynamic relation between the canonical and grand canonical ensemble. 
The canonical partition function with the system size N is given by 

. N N-l 

z N = ] n dd y l II f(v* - y i+1 ) = vtfm"- 1 (2-i) 

i=l i=l 

where V is the total volume of the system, and f(x) is a function assigned to each bond and 
damps sufficiently fast at long distance compared to the typical length scale a . f{p) is its 
Fourier transform: 

f(p) = fd d xe^f(x). (2.2) 



For example, we can take f(x) = exp(—(x/ao) 2 /2). 

Correlation functions for density p(x) = X^Li S^(x — y l ) can be easily calculated. One 
point function becomes 



1 N N-l n N 

Two-point function is defined as 

N N-l N N 

U dd y l II /(v i -v i+1 )E* (, V-i/ < )E 

J N J i=1 i=1 i= i j = i 

and its Fourier transformation is given by 
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where 

h(p) = f(p)/f(0)=exp(-(a p) 2 /2), (2.6) 

H (p) = l-h(p) = (a p) 2 /2 + --: (2.7) 

We can approximate the two-point function as 

9 ( n\p) = ^-777-n - T^T— + 0( f ~ vv> )) (2.i 



in the following scaling region: 
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(2.9) 



or 

N-^/ao < p < l/o . (2.10) 

The scaling region is between the ultraviolet cutoff scale a Q and the infrared scale of the 
system extent, which is given by £ = aoN 1 ^ 2 . 

The two-point correlation function is also calculated in the grand canonical ensemble. A 
conventional grand canonical partition function is given as 

Z K0 = £ k»Z n = -p^- (2.11) 

where Kq is fugacity and k c = (/(0)) _1 . For later convenience, we first define an unnormalized 
two-point function by 

oo 
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Hence the normalized two-point function becomes 



(2.12) 
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V(l-K f(0))(l-K f(p)) Vk c -k 
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where 8k = k c — kq. Since (N) = k c /Sk, the correlation function behaves as that of massive 
scalar particles 

»"<>>« ~ T aw + 1/ W (2 ' 14) 



and agrees with the result in the canonical ensemble calculation. This result gives the 
correlation length £ = aoN 1 ^ 2 , which indicates the Hausdorff dimension of random walk 
d H = 2. 

Here we give two different definitions of Hausdorff dimensions. The first one is defined in 
terms of the relation between the system size N and the extent of the system. The infrared 
behavior of the above two-point function shows that the correlation damps rapidly over the 
length scale £ = clqN 1 / 2 . Since the extent of the system L ~ £ is proportional to iV 1//2 , the 
Hausdorff dimension is given by d H = 2 where d H is defines as 

L = a N 1/d ». (2.15) 

The second definition is to use the behavior of the correlation function at much shorter 
length scale than the system size L — £. In rf-dimensional coordinate space, the density 
correlation behaves as 

gV(x) ~ ^ ( " m 2 N) (2.16) 

where m = l/£. If \x\ <C l/m, the mass term can be neglected and g^ 2 \x) ~ , } d _ 2 . The 
total number of points within a ball of radius R (R <^ L) around a certain point is evaluated 

as 

N(R) = J R d d xg^ {x)~(—\ (2.17) 

(2) (2) 

and gives the Hausdorff dimension d H = 2. Here the definition of d H is 

j( 2 ) 

/ R \ d n 

N(R) = (-) . (2.18) 

The second definition of the Hausdorff dimension is determined only through the behavior 
of correlation functions in the scaling region Oo < i < ^ and has nothing to do with 
the behaviors near the infrared cutoff. Hence, it is a more appropriate definition than the 
first one from the thermodynamic viewpoint. Of course, it is quite natural that we should 
obtain the same dimension d# from the both definitions since the extent of the system is 
approximately evaluated at the length L where N(L) = N. 

Finally in this section, we comment on generalized types of grand canonical ensembles. 
Since the motivation of introducing grand canonical ensembles is to reproduce the same 
thermodynamic quantities as those in the canonical one, we may assign different weights in 
summing over different N. We define generalized grand canonical ensembles by 



Z K(h i — 2^ N k q Z n , 

N=l 

oo 

G^Up) = E^C^(p), (2.19) 



N=l 



where kq is fugacity. The criterion for a 'good' grand canonical ensemble is such that we can 
take the correct thermodynamic limit, or, in other words, we can correctly take the universal 
part in the sum over N. That is, the correlation functions in the grand canonical ensemble 



at the critical value of fugacity should reproduce those in the canonical ensemble for large 
iV: 

lim^ 2) = lim s«. (2.20) 

AT— >oo «0— *«0,c u ' 

This criterion does hold if the grand canonical ensemble is dominated by large N systems. 
If the above criterion is satisfied for some value of I, it does hold for larger values of I. In 
the case of random walks, it already holds for the conventional grand canonical ensemble of 
I = and we do not need to introduce the generalized ensembles of I > 0. 

From eqs.( |2.1lD , ( |2.12| ), the generalized partition functions and the correlation functions 
are given by 

Z KOtl = ^ —j Z K0 , 

Op)=(^)^). (2.21) 

We are interested in the behaviors near the critical point 5k ~ (N —>■ oo) and in the scaling 
region where 1/(k c — k ) > 1/(k c — h(p)K ) ^> 1/k c . In this region, they become 
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+ ■■■). (2.22) 



Since (N) = Z KQi i + i/Z KOt i = k c (1 + 1)/Sk for the generalized grand canonical ensembles of I, 
the normalized correlation functions become 



which agrees with the previous result ( |2.14|) of / = 



3 Branched polymer dynamics 

3.1 Canonical Ensemble 

Branched polymers are a statistical system of N points connected by N — 1 bonds whose 
lengths are of order clq. The canonical partition function is defined as 



1 r N 

v, E /n<*V II fW-v'l (3-iJ 
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G:tree graph i=l (ij):bonf of G 

where f(x) is a function assigned to each bond in each graph, and it damps sufficiently fast 
at long distances compared to the typical length scale do- The presence of the factor 1/iV! 
is due to the fact that the N points are regarded identical. 



We can calculate a partition function for each N, by counting the number of all possible 
tree graphs: 
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z 2 -- 


= v, 


Z-i -- 
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In = ^N N - 2 f(0f- l V, (3.2) 



where / (p) is a Fourier transform of / (x) ; 

f(p)= [d*xe*°f(x), (3.3) 



and V is the total volume of the system. A derivation of the general form (|3.2j ) is given in 
the appendix A. 

We define an (unnormalized) m-point correlation function of density operators as 

^N \ x ) ■ ■ ■ 5 x ) 

= Z N < p(x 1 ) ■ ■ ■ p{x m ) > N 



1 N 

^ E /lI rf V n f(y l -v 3 ) p{x x )--- P {x m ), (3.4) 



G:tree graph i=l (ij):bond of G 

where the density operator is defined by 

N 



p(s) =£ <$«>(* -y<). (3.5) 

Due to the translational invariance, one point function is proportional to the partition func- 
tion; 

G® = yZ N . (3.6) 

The one-point function is nothing but the partition function with one marked point. 

3.2 Grand canonical ensemble 

We then define partition functions and m-point correlation functions in the generalized grand 
canonical ensembles as in section 2: 



Z K0 ,i = Y^N 1 k%Z Nj (3.7) 

N=l 

G%l(x\...,x m ) = J2N l K%GW(x\...,x™). (3.8) 



k is the fugacity. 

The criterion for a 'good' grand canonical ensemble is such that we can take the correct 
thermodynamic limit in the following sense. The correlation functions in the grand canonical 
ensemble at the critical value of fugacity should reproduce those in the canonical ensemble 
for large N: 

where we have defined normalized correlation functions as 

Mm) fi(m) 

(m) _ ^N (m) _ ^kq,1 /o i n \ 

9N — ~^ J 9 K0 ,l ~ ~~^ • {d.LV) 

In the case of random walks, we have confirmed that it does hold for any nonnegative value 
of I but we need to check it in the case of branched polymers. To satisfy this criterion, 
the grand canonical ensembles with iV-dependent weights should be dominated by large N 
systems. This is not assured only by taking the fugacity near the critical value. This is 
because, if G™ behaves as (k 0iC )~ n N a for large N and I + a < — 1, the summation over 
iV is dominated by small N system, not by the large N ~ k c /Ak even near the critical 
point. On the other hand, if we take a sufficiently large /, we can obtain the correct large iV 
correlation functions in the grand canonical ensembles, which are, of course, independent of 
I. 

We illustrate the above mentioned criterion by taking the partition function as an exam- 
ple. Since the canonical ensemble partition function ( |3.2|) behaves at large N as 

AT-5/2 

Z N ' 



'2ne- N 

the grand canonical ensemble is approximated by 

V °° 
Z KO i ~ ■==* ViV*- 5 / 2 (k/k c 
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where 



k = /(0)Ko, 



Kr. — 6 



An = k c — k. (3.13) 



If we take I sufficiently large, the integrand in eq. ( [3.12 ) has a peak at N ~ k c /Ak and we 



can make N large by letting k approach k c . On the other hand, if / is not sufficiently large, a 
non-universal small iV behavior dominates the summation and we cannot obtain the correct 
answer of the large iV limit by a grand canonical ensemble. 



3.3 Schwinger-Dyson eq. 

In this subsection, we recapitulate the arguments that the correlation functions for branched 
polymers in the conventional grand canonical ensemble are given by massless 3 theory. 
Let's consider the following correlation functions G^ = G™ 1=0 , which are suitable for 
Schwinger-Dyson analysis: 



/~t(yn) i ~X 
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We write a Fourier transform of G^ (x 1 
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Schwinger Dyson equation for 1-point function GW becomes 



where 



h = Ke h 



b = f(0)&£, 



as can be seen from figure [I]. Figure g illustrates eq. ( ft.lfij ) 



(3.14) 



(3.15) 

(3.16) 
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Figure 1: Schwinger-Dyson equation for one point function. A grey blob and a black point 
mean G^ and Kq, respectively. 



At the critical point, 



k c — e , 
db/dk diverges. Near this critical point, N becomes large; 

Ab ~ \F2e\TK^ ~ 1/VN, 



(3.18) 



(3.19) 




Figure 2: Schwinger-Dyson equation, k = be b . At the critical point, b c = 1, k c = e x . 



where 



Ab = b c - b, 



(3.20) 



The one-point function (which is equal to the partition function of I = 1) now behaves as 
follows. 



Z Kn i 



K0,l = l 



N N 

bV V 



/(0) /(O) 



1 - V2eVA«). 



(3.21) 



Next, we consider the 2-point function G^(p). When we pick up any two points on a 
tree graph, we can fix the path connecting these two points. Thus, as can be seen from figure 
KL 2-point function is calculated to be 



Gff(p) = T,f(pnGg 



)Y+i 



s=0 



/(0)(1-W»(p)) 



(3.22) 



where 



h(p) = /(p)//(0) = 1 - H{p) 



2„2 



1 — ca p 



(3.23) 



Here c is a positive constant of order one. Recall that f{x) damps rapidly out of the region 
< x < do- 
Near the critical point, b ~ b c = 1, the 2-point correlation function behaves as 



g£(p) 



i 



/(0) H(p) + Abh(p) 

1 
H(p) + N-V 2 ' 



(3.24) 



A 

G ,2 = 




Figure 3: 2-point function G^ is made out of 1-point function G)£f, which is written by a 
gray blob. 



Here we used eq. ( |3.i9| ). Thus, the correlation length is £ = a^N 1 ^, which shows the first 
definition of the Hausdorff dimension defined in (2.15) to be d H = 4. Let us consider the 
following region: 

a < x < £ = a N 1/4 (3.25) 

or in momentum space 

I = N'^/ao < p < l/o . (3.26) 

ao gives an ultraviolet cut-off whereas £ gives an infrared cut-off length over which correlation 
functions damp rapidly. In the above scaling region of Q3.26| ), the correlation function behaves 



as an ordinary massless field G^(p) ~ 1/p 2 and gives the second definition of the Hausdorff 
dimension of ( 2.1 8|) , d H = 2. This is different from the Hausdorff dimension d H determined 



^from the relation of the system size N and the extent of the system. In the next section, we 
show that the above behavior of the correlation function in the scaling region is not correct 
and hence gives the incorrect Hausdorff dimension d H . 

Finally, we consider m > 2 point correlation functions. As in the case of the two-point 
function, when m points are fixed on each tree graph, we can uniquely fix the path connecting 
these points. Therefore, an m-point function &™> is represented as a summation over all 
tree diagrams with m fixed points in which G^J (p) appear as propagators. For example, 



Gg(p,g) = (GU)- 2 GW(p)G$( q )G$(p + q ). (3.27) 

In general, we obtain the following result for m-point correlation functions; 

G^™' ~ correlation functions of massless 3 theory at tree level. (3.28) 

In the next subsection, we consider generalized grand canonical ensembles with larger I and 
we point out that the results ( |3.28|) do not correspond to the correct thermodynamic limit. 



3.4 Correlation functions in thermodynamic limit 

Let us consider the generalized m-point correlation functions with I > 1. From the definition 
T8j) , they can be obtained by applying Z-th derivative to the I = case: 

GfflCpV--^- 1 ) = (ko^^^Cp 1 ,...,^ 1 ) 



) l GtU{p\---,P m - X ) (3-29) 



l-bdb' K0 ' l=0 
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Partition functions are obtained from fl3.21| ) as 



7 - ( b g V-i W 

^ ,i y 1 _ hdh > ^ (0) 

Near the critical point (b c — 1), they become 

V (2/ -5)! 
/(0)(l-6) 

2-point function with / = 1 is given by 



(3.30) 



Z «o,«>2 ~ TTTTT^ E^T¥- ( 3 - 31 ^ 
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GZUv) = (t—^)( 



i-bdb jy f(o)(i-bh( P )y 
i i 

/(0)(l-6)(l-&%)) 2 

1 

~4 



p 



(3.32) 



It is because in the scaling region ( |3.26| ) and near the critical point, the following inequality 
holds: 

b~ l<< T^bh(f) ~ c{a, P f + iV-V 2 « 1^6 ~ ^ ( 3 - 33 ) 

This behavior is different from that of G K( j l=0 ~ 1/p 2 . Similarly for / > 2, the behavior of 
2-point function becomes 



, (2 ) J. (21 -3)!! 1-6 1-2 1-6 . 

<Ko,i>2(P)~ /( )(i-6)a-i (l-&%))2 [ + 1-6%) + 2l-3 ( l-i%) 



Due to the inequality ( |3.33|) , the derivative J^ acts dominantly on j^, not on 1 _ fe 1 fe( r . The 



normalized 2-point functions now become 

a(2 ) n 1 21-3 1 Q (l-fr) 2 1 v 

9ko,i>2{P) ~ y(l-6) 2 (l-%))2 ( 2Z-3 (l-%)) 2 ' j 

~ ^^ (1 - 3 m^ +0(( w^ )2)) ' (3 ' 35) 

Here we have used (iV) = Z KOi i + i/Z Ko j = (2/ — 3)/(l — b) 2 for the generalized grand canonical 
ensembles of /. Their p-dependences are all the same except / = case. The I = case, which 
can be obtained directly from the Schwinger- Dyson equation, does not reproduce the correct 
thermodynamic result. Instead, we should consider a 'good' grand canonical correlation 
function with / > 1, otherwise a non-universal small N behavior affects the summation 
and we cannot obtain the universal result. In the appendix B, we estimate the large N 
asymptotic behavior of the two-point function by the saddle point method. Such an explicit 
analytical result is completely consistent with the analysis here. 

The behavior of g L>1 (p) ~ 1/p 4 gives (the second definition of ) the Hausdorff dimension 



( 2.18| ) d H = 4, which is now consistent with (the first definition of) the Hausdorff dimension 

11 



discussed in the previous subsection. An argument expected from the figure [3| is that the 
effect of branching could be absorbed by renormalizing the mass. If so, the propagator would 
behave as that of random walks with a renormalized mass and we might obtain the identical 
result with the / = case. We discuss in the next section why this argument is not correct. 
Similarly, 3-point functions become 

9k1i=o(p^) ~ g{p)g{q)g{p + q), 

<S>i(p, q) ~ g(p)'g(q)g(p + q) + g(p)g(i) 'g(p + q) + g(p)g(q)g(p + q)', (3.36) 



where 



g{p) 



1 — bh(j>) p 2 



g ' ip) = (1 - bh(p)Y ~ J' (3 ' 37) 

The behaviors do not change above 1 = 1. Only one propagator in a graph is replaced by 
g'(p). This is because the derivative Jg dominantly act on the factor 1/(1 — b) rather than 
on 1/(1 — bh(p)), as in the case of 2-point functions. Here again, we might be tempted to 
argue that the only effect of branching is mass renormalization and to give rise to a three 
point vertex. If so, we might conclude that the 3-point function behaves as the result of 
/ = case. As we discuss in detail in the next section, this argument is not correct in the 
thermodynamic limit and we should only retain I > 1 cases. Hence correlation functions 
for branched polymers are expressed in terms of 3 theory at tree level with a single mass 
insertion. 

For m-point functions (m > 3), we can obtain the same result. 



y^i = o ~ correlation functions for (ft 3 theory at tree level 

gi^i>i ~ correlation functions for 3 theory at tree level with a mass insertion. 



(3.38) 

The universal correlation functions with I > 1 represent the correct correlation functions in 
the thermodynamic limit. 

As a consistency check, the following relation between an (m+l)-point function and an 
m-point function must hold: 

&i(p\ ■ ■ -,P m -\p m = 0) = Ng^(p\ ■ ■ • ,p- 1 ). (3.39) 



It actually holds because in the L.H.S. of eq. (|3.39|) , the special class of diagrams dominate 



in which the m-th end point is attached to a propagator g'(p). It is due to the inequality 



g'(p — fyg(p) ^ g'(p)g(p — 0)- Therefore, it is equal to the R.H.S. of eq. (|3.39| ) 



4 Physical interpretation by a single mother universe 

In this section, we give a physical interpretation why the propagator behaves as 1/p 4 instead 
of 1/p 2 . As we can see from the figure |3|, the effect of branching seems to be absorbed by 
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renormalizing the mass and we might conclude that the two-point function behaves as that 
of random walks. If this is the case, the propagator should be given by an ordinary massive 
scalar particle with a renormalized mass (ao-/V 1//4 ) _1 , instead of (ao-ZV 1 / 2 ) -1 . Similarly, higher 
point correlation functions should be given by tree graphs of 3 field theory. As we saw in 
the previous section, these are not the correct behaviors of correlations. In this section, we 
give a physical interpretation why the propagator behaves as 1/p 4 instead of a conventional 
behavior 1/p 2 and why the higher point correlation functions behave as 3 theory with a 
single mass insertion. 

Let us consider the two-point function as an example. Two-point function G K l (p) was 
defined by a sum of two-point functions in the canonical ensemble (eq. (|3.14|) ) : 



G ( S(p)=J:^gP( P ). (4.1) 



MM 2 ) 

N=l 



On the other hand, from the equation ( |3.22| ), it is written as a sum of all contributions over 
s where s is the length between the two points in concern. Each term G^ represents a 
gray blob in figure |3|. When we fix the total number of the system N, the N points are 
distributed among (s + 1) blobs. We will show here that the most dominant contributions to 
the correlation functions are those that most of TV points are concentrated on only a single 
blob. We call this blob the mother universe. In the branched polymer, there is only one 
mother universe and the other universes (blobs) contain much fewer points than the mother 
universe. 

To show this, we first note that the 1-point function (blob) is expanded as 



GS = £ ^Un (4.2) 

n=l 



where u n is given at large n as 



„„ = G m „ "-y (4.3) 



By using this expansion and eq. (|3.22|) , we obtain the two-point function for fixed N as 



oo oo 

.(2), 



G { n\p) = £ f(PT £ • • • £ «no • • • U n Mno + -n s ■ (4.4) 

s=0 yno=l n a =l J 

Each contribution in the bracket comes from a graph in which the first blob contains hq 
points, the second n\ points, and so on. Using ( |4.3| ), the term in the bracket becomes 

4= /(0)"— V £ • - £ n-^W ■ ■ n^5 N , no+ ... ns . (4.5) 

V v 2tt J 



no=l n a =X 

In the case of s — 1, the summation 

N-l 

Y^n-^iN-n)-^ 2 . (4.6) 



n=l 
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is dominated by terms of n ~ and n ~ N '. For the general case with the exponent a, 

N-l 

Y, n a (N - n) a , (4.7) 

n=l 

the sum is dominated at the boundaries for a < — 1 and asymmetry arises between two 
blobs. On the contrary, for a > — 1, N points are distributed equally and neither blob is 
special. This argument can be generalized to s > 1. To conclude, most of N points belong to 
a single blob along the propagator. We then have to divide the propagator with length s at 
the mother universe. Since the other blobs contain only a finite number of points, the effect 
of branching other than dividing the propagator into two pieces is simply to renormalize the 
mass of propagator. Hence, the propagator behaves as a product of two ordinary ones with 
a renormalized mass; 

\p z + m A J 

We can also apply a similar argument for the higher-point functions. If the effect of 
branching is only to renormalize the mass term, the higher-point correlation functions will be 
represented by diagrams of 3 theory with propagators l/(p 2 + m 2 ), which is not the correct 
thermodynamic behavior of the correlation functions as we saw in the previous section. 
Similar to the case of 2-point functions above, we can argue that there is only one mother 
universe in which most of the N points reside. Since only a single blob contains an infinitely 
many points, we have to divide one of the propagators at the mother universe and this 
propagator behaves as l/(p 2 + m 2 ) 2 . This mother universe corresponds to the mass insertion. 
The other blobs contain a finite number of points and the effect of branching can be absorbed 
into mass renormalization. This is the physical reason why the higher point correlation 
functions are represented by tree diagrams of 3 field theory with a single mass insertion. 

The above statement that there is only one mother universe in the branched polymer is 
explained differently as follows. Let us consider again the original branched polymer systems 
with iV points and N — 1 bonds. By counting the number of ways in which we can divide a 
branched polymer into two parts by cutting a bond, we obtain a relation 

N-l 

{N-l)Z N /V = f{0)Yl [N'Z N ,/V\ [{N-N')Z {N „ W) /V\. (4.9) 

JV'=1 

The factor (N — 1) in the L.H.S. is interpreted as the number of bonds we can cut to divide 
the whole into two parts. The factors N' and (N — N') in the R.H.S. are interpreted as the 
number of points to which the bond connecting the two parts is attached. Since Z^ behaves 
as in eq.( |3.11| ) at large N, 

iV(l-e) 

NZ N /V^>f(0) J2 WZ N ,/V} [{N - N')Z {N -N')/Vi (4.10) 

N'=Ne 



which means the summation in the R.H.S. of eq.( O ) is dominated by terms of N' ~ or 



N' ~ N. This formula can be interpreted as follows. If we divide any graph into two parts 
by cutting some bond, we find only finite points in one of them and most of them belong 
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to the other. The mother universe belongs to the larger part. Consider the 2-point function 
as an example. By dividing the graph of figure |] into two parts at some bond, we find that 
most points are dominantly distributed on only one of them. 

We can apply the same procedure to the dominant part with infinite points. After 
repeating it several times to divide the total graph into several pieces, we find that only a 
single part consists of infinitely many points and the others consist of finite points. We can 
find out which blob in figure || is the mother universe when, after the several repetitions, the 
dominant part with infinite points is detached from the path in figure |3|. We can also apply 
the same argument to the higher-point functions. Here we point out the difference of the 
argument given here and that for random surfaces ||. Although a similar inequality appears 
for random surfaces, the factor N is absent on the L.H.S. of the inequality Q4.10D in that 
case. 

5 Conclusion and discussion 

In this paper we have shown that the correlation functions for branched polymers are given 
by those for 3 theory at tree level with a single mass insertion if we correctly take the 
thermodynamic limit. It is not given by those for 3 theory at tree level themselves. We have 
interpreted the single mass insertion as the presence of a single mother universe in branched 
polymers. 

We have reviewed random walks in order to clarify the relation between the canonical 
and grand canonical ensemble. We have introduced generalized grand canonical ensembles 
which assign different weights for different system sizes. We have emphasized that a 'good' 
grand canonical ensemble is such that the ensemble average should be dominated by the 
systems of large size. In branched polymers, this criterion is not satisfied in the conventional 
grand canonical ensemble. Nevertheless we can consider 'good' grand canonical ensembles 
in branched polymers. Our conclusion follows as the universal prediction of 'good' grand 
canonical ensembles. It represents the correct scaling behavior of the correlation functions 
in canonical ensembles of large system size N. 

Acknowledgments 

This work is supported in part by the Grant-in- Aid for Scientific Research from the Ministry 
of Education, Science and Culture of Japan and by the National Science Foundation under 
Grant No. PHYS94-07194. 

Appendix A 

In this appendix, we derive the canonical ensemble partition function (|3.2|) from the Schwinger- 
Dyson equation ( |3.1(j| ). Let us solve b as a form of expansion in k. Each coefficient is 
calculated to be 

; dn^^ —I db(l-b)e~ b Tl — t-^-t (A.l) 

(A.2) 
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2iriJ K =o k n+1 2iriJb=o ' (be~ b ) 



Hence, 

OO ATM — 1 

jv=i iV - 



From eqs. ( |3.6|) , (|3.14|) and (|3.17| ), b is expanded as 



00 N 
b = /(0) E 17 ^ 4- (A.4) 

iV=l K 



Comparing these two expansions, we get the result of eq. 

Appendix B 

In this appendix, we derive the two-point correlation function in the canonical ensemble for 
fixed but large N. Similar to the calculation in appendix A, we can obtain the two-point 
function for a fixed N from eq. Q3.22j ); 



lltl L=0 K N+1 

1 r „ 1 1 - b 



G ( n\p) = ^~.j>dn 



J(0)J 27TZ Jk=0 k 



1 db^ — — -b~ N e bN . (B.l) 

2m Ib=o /(0) 1 - bh{p) 

This integration can be estimated for large N by the saddle point approximation. Since 
b~ N e bN = e N ( b - l °g b )^ the saddle point is at b = 1. The steepest descent direction is along the 
imaginary direction. We change the variable from b to t around the saddle point as 

b-\ogb=l-^. (B.2) 

The new parameter t is written in terms of b by 

t = z(l-6)(l + i^ + ^(l-6) 2 + ---). (B.3) 

Solving this, we obtain 

1 9 13 

--i = -*t(i--*t-^; 2 + ■•■)• (B.4) 

Therefore in the saddle point approximation, the correlation function becomes 

[f(0)J N [P) ~ 2m J-Jl/b-l + Hipf 

e^ r ,, t 2 (l-13t 2 /36 + ---) _ m , /2 

2vri-oo (H(p)-2t 2 /3 + ---) 2 + t 2 (l-13t 2 /36 + ---) 2 



2tt J -00 (H(p) - 2t 2 /3N + • • -) 2 + t 2 /N(l - 13t 2 /3GN + • • -) 2 
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where H(p) = 1 — h(p) ~ c(aop) 2 . Hence the normalized correlation function becomes 
^(P)- 6S>(P) 



Z 



N r°° _dt t 2 (l-13t 2 /36iV + ---) _ f2/2 

V 7-00 v 7 ^ (#(p) - 2t 2 /3N + • • -) 2 + t 2 /AT(l - 13t 2 /36iV + •••)** { ' ' 

H [°° _^_ t p -tV2 M , oi lx, m 7) 

V J-oo ^H(p) 2 -U 2 H(p)/3N + t 2 /N l k A^' l ' ; 



The correlation function for branched polymers behaves very differently from the case of 
simpler random walks. In the scaling region 1/y/N < H(p) <^i 1, the 2-point function can 
be evaluated as 

*(2), s N f°° dt t 2 .2 in .„ , 

9 " (p) ~ vLvS hWTWJn * ' < R8) 

This result for canonical ensemble is completely consistent with the result of the grand 
canonical ensemble Q3.35p . 
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